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Synopsis 

In plants, sulfur must be obtained from the environment and assimilated into usable forms for metabolism. ATP 
sulfurylase catalyses the thermodynamically unfavourable formation of a mixed phosphosulfate anhydride in APS 
(adenosine 5'-phosphosulfate) from ATP and sulfate as the first committed step of sulfur assimilation in plants. In 
contrast to the multi-functional, allosterically regulated ATP sulfurylases from bacteria, fungi and mammals, the plant 
enzyme functions as a mono-functional, non-allosteric homodimer. Owing to these differences, here we examine the 
kinetic mechanism of soybean ATP sulfurylase [GmATPSl {Glycine max (soybean) ATP sulfurylase isoform 1)]. For 
the forward reaction (APS synthesis), initial velocity methods indicate a single-displacement mechanism. Dead-end 
inhibition studies with chlorate showed competitive inhibition versus sulfate and non-competitive inhibition versus 
APS. Initial velocity studies of the reverse reaction (ATP synthesis) demonstrate a sequential mechanism with global 
fitting analysis suggesting an ordered binding of substrates. ITC (isothermal titration calorimetry) showed tight binding 
of APS to GmATPSl. In contrast, binding of PPj (pyrophosphate) to GmATPSl was not detected, although titration of 
the E»APS complex with PPj in the absence of magnesium displayed ternary complex formation. These results 
suggest a kinetic mechanism in which ATP and APS are the first substrates bound in the forward and reverse 
reactions, respectively. 
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INTRODUCTION 



The most abundant environmental source of sulfur is sulfate 
(S0 4 2_ ), which is a chemically inert molecule [1,2]. For plants 
and microbes to utilize this essential nutrient, sulfate is enzymat- 
ically converted into a chemical species that is energetically 
favourable for reduction [3,4]. The sulfur assimilation pathway 
provides sulfide for a range of biosynthetic pathways that 
supply methionine, glutathione, iron-sulfur clusters, vitamin 
cofactors such as biotin and thiamin, and multiple specialized 
metabolites such as glucosinolates [1,2,5,6]. 

The first enzymatic reaction in the sulfur assimilation pathway 
of plants is the non-reductive adenylation of sulfate catalysed by 
ATP sulfurylase (ATP: sulfate adenylyl transferase; EC 2.7.7.4) 
to yield APS (adenosine S'-phosphosulfate) and PPi (pyrophos- 
phate) (Figure 1) [7]. Generation of the mixed phosphosulfate 
anhydride bond in APS yields a high-energy molecule that drives 



subsequent reduction reactions in the assimilatory pathway [8- 
11]. APS formation is energetically unfavourable (^f e q~10 -7 - 
10 ~ 9 ) compared with the reverse reaction and requires high sub- 
strate concentrations and downstream enzymes to maintain the 
forward reaction equilibrium in the pathway [7,8]. Downstream 
enzymes in the plant sulfur assimilation pathway (i.e. APS re- 
ductase and sulfite reductase) convert APS into sulfite and then 
sulfide for subsequent incorporation into cysteine [1-4,6,12]. 
Alternatively, phosphorylation of APS to PAPS (adenosine 3'- 
phosphate S'-phosphosulfate) supplies a sulfur-donor molecule 
for a variety of biosynthetic pathways [13-15]. In addition, coup- 
ling of the ATP sulfurylase reaction with hydrolysis of PPi by 
pyrophosphatase helps to maintain APS synthesis in vivo [16]. 

In plants, gene families encode multiple isoforms of ATP 
sulfurylase with varied expression patterns and organelle loc- 
alization [4,10,17-23]. The different genes encode plastidic and 
cytosolic forms of the enzyme in Arabidopsis thaliana (thale 
cress), Brassica juncea (Indian mustard), Solanum tuberosum 



Abbreviations used: APS, adenosine S'-phosphosulfate; GmATPSl, Glycine max (soybean) ATP sulfurylase isoform 1; ITC, isothermal titration calorimetry; PAPS, adenosine 
S'-phosphate S'-phosphosulfate; PPi, pyrophosphate. 

1 To whom correspondence should be addressed (emailjjez@biology2.wustl.edu). 



© 2013 The Author(s) This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC-BY) (http://creativecommons.Org/licenses/by/3.0/) 585 
which permits unrestricted use, distribution and reproduction in any medium, provided the original work is properly cited. 



G. E. Ravilious and others 



04,0- "CX« Q,* 

sulfate f, pyrophosphate ° " 

O 0 0 

°o rv 

+ o + °< 

ATP > C 0 APS O, ( 

XjTZ Sulfurylase hJ, ^ 

^ HO OH 

Figure 1 Overall reaction catalysed by ATP sulfurylase 



(potato) and Glycine max (soybean) and share >60% amino 
acid sequence identity [4,17,18,20,22]. For example, the soybean 
genome contains four ATP sulfurylase genes (Glymal0g38760, 
Glyma20g28980, Glymal3g06940 and Glymal9g05020) encod- 
ing proteins with predicted plastidial and/or mitochondrial local- 
ization [4]. 

Biochemical analysis of a plastidic ATP sulfurylase from soy- 
bean [GmATPSl {Glycine max (soybean) ATP sulfurylase iso- 
form 1)] showed that the enzyme functions as a homodimer [10]. 
The oligomeric structure of the plant ATP sulfurylase differs from 
that of the enzyme from bacteria, fungi and mammals and also 
lacks allosteric regulation. Prokaryotic forms of ATP sulfurylase 
are heterodimeric proteins in which a GTPase subunit allosteric- 
ally activates the catalytic subunit [24]. In fungi, such as Sacchar- 
omyces cerevisiae and Pencillium chrysogenum, ATP sulfurylase 
functions as a homohexamer in which each monomer contains an 
N-terminal ATP sulfurylase domain and a regulatory C-terminal 
APS kinase domain [25,26]. The human enzyme functions as a 
homodimeric protein with a reversed ATP sulfurylase/ APS kinase 
domain order compared with the yeast enzyme [27]. In fungi and 
mammals, the fusion of two activities into a single protein al- 
lows for direct synthesis of PAPS, which is used as a substrate 
for subsequent sulfonation reactions [28]. Interestingly, in certain 
diatom and microalgae species, ATP sulfurylase and pyrophos- 
phatase are also reported to function as a fusion protein [16]. 
In addition, potential interaction between ATP sulfurylase and 
APR reductase as a protein complex in onion suggests additional 
organization of enzyme activities in the plant sulfur assimilation 
pathway [29]. 

Owing to the differences in oligomeric organization, domain 
architecture and biochemical regulation between the ATP sul- 
furylases from plants and other organisms, here we investigate 
the kinetic mechanism of a plant ATP sulfurylase (i.e. GmATPS 1 ) 
to understand the reaction sequence catalysed by this critical en- 
zyme of the sulfur assimilation pathway in plants. 



EXPERIMENTAL 

Protein expression and purification 

The pET-28a-GmATPSA48 bacterial expression construct, 
which encodes GmATPSl as an N-terminal hexahistidine-tagged 
protein lacking the plastid localization sequence, was previously 



described [10]. Transformed Escherichia coli BL21(DE3) were 
grown at 37 °C in Terrific broth containing 50 /xg ml ~ 1 kanamy- 
cin until A 60 o n m~ 0.6-0.9. Following induction with 1 mM isop- 
ropyl 1-thio-yS-D-galactopyranoside, cell cultures were incubated 
overnight (20 °C; 250 rev./min). Cell pellets were prepared by 
centrifugation (10 000 g; 15 min) and then resuspended in 50 mM 
Tris (pH 8.0), 500 mM NaCl, 20 mM imidazole, 10% (v/v) gly- 
cerol and 1 % (v/v) Tween-20. Cells were ruptured by sonica- 
tion and the lysate cleared by centrifugation (45 000 g; 30 min). 
The lysate was loaded onto an Ni 2 + -nitriloacetic acid (Qiagen) 
column. After extensive washing with lysate buffer minus Tween- 
20, the bound protein eluted using 250 mM imidazole in wash 
buffer. Incubation with thrombin during overnight dialysis at 4 °C 
against wash buffer removed the His-tag. Dialysed protein was re- 
loaded on a mixed benzamidine-Sepharose/Ni 2 + -NTA column. 
The flow-through was loaded onto a Superdex-200 26/60 Hi- 
Load FPLC size-exclusion column equilibrated in 25 mM Hepes 
(pH 7.5), 200 mM KC1, and 5 % glycerol. Fractions correspond- 
ing to the major protein peak were pooled and examined for purity 
by SDS-PAGE. Protein was concentrated (Amicon) to 5-10 mg 
ml ~ 1 with protein concentration determined using a molar ex- 
tinction coefficient (£280nm = 54890 M~ 1 cm -1 ) calculated with 
ProtParam (http://web.expasy.org/protparam). Protein was flash 
frozen in liquid nitrogen and stored at — 80 °C until used. 

For assays of ATP sulfurylase activity in the forward (i.e. APS 
synthesis) direction, A. thaliana APS kinase was used. Details 
of the generation of the expression construct, E. coli expression, 
protein purification, and enzyme assay for AtAPSK (A. thaliana 
APSK) were previously described [10,30-32]. 

Enzyme assays 

Initial reaction velocities were determined by observing the 
rate of change in absorbance of pyridine nucleotide at 340 nm 
(e = 6270 M _1 cm -1 ) in 500 \x\ systems at 25 °C using a Beck- 
man DU800 UV/vis spectrophotometer. The forward APS syn- 
thesis reaction used an assay system consisting of 50 mM Tris 
(pH8.0), 15 mM MgCl 2 , 100 mM NaCl, 0.4 mM phosphoen- 
olpyruvate, 0.2 mM NADH, 0.05 units of APS kinase, 20 units 
of pyruvate kinase and 30 units of lactate dehydrogenase. The 
specific activity of APS kinase (1.5/xmol min -1 mg -1 ) was 
determined spectrophotometrically, as described elsewhere [10]. 
The reverse ATP synthesis reaction used an assay system of 
50 mM Tris (pH 8.0), 5 mM MgCl 2 , 1 mM NADP+ , 1 mM gluc- 
ose, 2 units of hexokinase and 1 unit of glucose-6-phosphate de- 
hydrogenase. All reactions were initiated by addition of enzyme 
and were corrected for non-enzymatic rates. 

For analysis of the bi bi substrate kinetic mechanism of 
GmATPSl, initial velocity rates were measured under standard 
assay conditions with a matrix of substrate concentrations. In 
the forward direction, assays used varied Na 2 S0 4 and ATP con- 
centrations. For the reverse reaction, initial velocities were de- 
termined using varied APS and sodium PPi concentrations. The 
resulting data were analysed by global curve fitting in SigmaPlot 
(Systat Software, Inc.) to model the kinetic data to rapid equi- 
librium rate equations describing ordered sequential, v = (V max 
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[A] [B])/(K A K B +K B [A] + [A] [B]), and random sequential, 
v = (V max [A] [B])/(a K A K B + K B [A] + K A [B] + [A] [B]), kin- 
etic mechanisms, where v is the initial velocity, V mSLX is the max- 
imum velocity, K A and K B are the K M values for substrates A 
and B, respectively, and a is the interaction factor if the binding 
of one substrate changes the dissociation constant for the other 
[33]. 

For the inhibition studies, initial velocities were determined in 
reactions containing chlorate (0-0.3 mM) and varied concentra- 
tions of either ATP or Na 2 S0 4 in the forward reaction. The global 
fitting analysis of inhibition experiments used SigmaPlot to fit all 
data to the equations for competitive, v = V max /(\ + ((K m /[SJ) 
(1 + [I\/K i} )) or uncompetitive, v = V m J(l + [1]/^ + KJ[S]) 
inhibition [34]. 



Calorimetric measurements 

ITC (isothermal titration calorimetry) experiments were per- 
formed using a VP-ITC calorimeter (Microcal, Inc.). GmATPSl 
was dialysed at 4°C in 25 mM HEPES, pH 7.5, 200 mM KC1, 
5 % glycerol. Stock solutions of APS and PPi were prepared in 
the same buffer. Protein and ligand solutions were degassed at 
room temperature (22 °C) prior to use. For each titration, 20-30 
injections of 10 /xl of ligand were added into sample solutions 
containing protein (18.5-50 /xM) in the presence or absence of 
ligands. GmATPSl complexed with APS was formed by incub- 
ating protein and 0.1 mM APS 30 min (4°C) followed by equi- 
libration at 25 °C before titration. Data were analysed using a 
one-site (i.e. identical sites) binding model, as follows: 

g; ot = Vo • Mj ot • ((nK^AH^/l + K x x), 

where g/ tot is the total heat after the /th injection, V 0 is the 
calorimetric cell volume, M,- tot is the concentration of protein in 
the cell after the /th injection, AH is the corresponding enthalpy 
change, n is the number of nucleotide-binding sites per monomer, 
and K is the equilibrium-binding constant. Fitting of data was 
performed with Origin software. 



RESULTS 



Forward APS synthesis reaction: initial velocity 
analysis 

The kinetic mechanism of the physiological APS synthesis re- 
action of GmATPSl was first examined using a bi substrate 
variation experiment. For these experiments, GmATPSl was ex- 
pressed in E. coli as an N-terminally hexa-histidine tagged protein 
and was purified by nickel- affinity and size-exclusion chromato- 
graphies [10]. Three possible kinetic models - ping-pong, ordered 
sequential and random sequential describe two substrate to two 
product or bi bi kinetic mechanisms [33]. 

Initial velocity data generated by co-varying the concentra- 
tions of ATP and Na 2 S0 4 was used to distinguish between the 
three possible mechanisms (Figure 2). The primary pattern, in 
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Figure 2 Initial velocity variation of substrates in the forward 
reaction 

Experimental data are shown as symbols in the double-reciprocal plots 
and lines indicate the global fit of all data to the equation for an ordered 
sequential bi bi mechanism. (A) Double-reciprocal plot of 1/v versus 
1/[ATP] at 0.5, 1.0, 2.5 and 5 mM sulfate (top to bottom). (B) Double-re- 
ciprocal plot of 1/v versus l/[sulfate] at 0.3, 0.5, 0.57 and 1.2 mM 
ATP (top to bottom). 



Table 1 Kinetic constants for an ordered sequential bi bi 
kinetic model 



Reaction 



Fitted parameters 



Forward 

Vmax (^mol min -1 mg of protein -1 ) 
K m ATP (fiM) 
K m sulfate (fiM) 
Reverse 

Vmax (Atmol min -1 mg of protein -1 ) 
K m APS (/xM) 

^pyrophosphate (/xM) 



14.6 ±3.5 
221 + 41 
1,030 ±280 

58.8 ±11.5 
28.3 ±0.8 
76.2 ±7.4 



which the lines intersect, indicates that GmATPSl uses a single 
displacement type of mechanism and eliminates the ping-pong 
type of mechanism from further consideration. The resulting 
data was globally fit to the equations describing the random 
and ordered bi bi substrate mechanisms. Examination of double- 
reciprocal plots, residual errors, standard errors of the fitted para- 
meters, and the correlation coefficient evaluated quality of the 
fit to the observed data. Modelling of the initial velocity data to 
a random sequential mechanism yields a correlation coefficient 
(r 2 = 0.989) slightly better than the fit to an ordered mechanism 
(r 2 = 0.980), but the interaction factor (a) for the random mech- 
anism was 9.53, which suggests that binding of one substrate is a 
negative interaction that decreases the affinity for the second sub- 
strate. Also, as discussed below, dead-end inhibition studies were 
consistent with an ordered mechanism. Thus, Table 1 summarizes 
the fitted kinetic parameters for an ordered sequential bi bi kin- 
etic mechanism. The fitted kinetic parameters reported here were 
comparable with those previously published for GmATPSl [10]. 



Forward APS synthesis reaction: dead-end inhibition 

To further examine the order of substrate addition in the forward 
reaction, chlorate (C10 3 ~) was used as a dead-end inhibitor. 
Chlorate has also been demonstrated in tissue culture to reduce 
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Figure 3 Chlorate inhibition in the forward reaction 

Experimental data are shown as symbols in the double-reciprocal plots 
and lines indicate the global fit of all data to the equation for either 
uncompetitive (A) or competitive (B) inhibition. (A) Double-reciprocal 
plot of 1/v versus 1/[ATP] at 0.3, 0.2, 0.1 and 0 mM chlorate (top 
to bottom). (B) Double-reciprocal plot of 1/v versus l/[sulfate] at 0.3, 
0.2, 0.1 and 0 mM chlorate (top to bottom). 



Figure 4 Initial velocity variation of substrates in the reverse 
reaction 

Experimental data are shown as symbols in the double-reciprocal plots 
and lines indicate the global fit of all data to the equation for an 
ordered sequential bi bi mechanism. (A) Double-reciprocal plot of 1/v 
versus 1/[APS] at 6.25, 12.5, 25 and 50 /xM PP, (top to bottom). (B) 
Double-reciprocal plot of 1/v versus 1/[PP\] at 1.25, 2.5, 3.0, 4.0, 5.0 
and 10.0 /xM APS (top to bottom). 



sulfur assimilation pathway activity [35] and is a dead-end inhib- 
itor of ATP sulfurylase from P. chrysogenum [36]. This approach 
was chosen instead of product inhibition studies because of the 
use of APS kinase in the coupling assay [8]. 

The effect of oxyanions on the forward reaction catalysed by 
GmATPSl was tested in the presence of chlorate, nitrate (N0 3 ~ ) 
and carbonate (C0 3 ~ ). Nitrate and carbonate were chosen be- 
cause they are naturally present in the chloroplast and may reg- 
ulate ATP sulfurylase activity [37]. In the presence of 100 /xM 
chlorate, nitrate or carbonate, GmATPSl displayed 35, 66 and 
80% of wild-type APS synthesis activity, which is consistent 
with oxyanion inhibition reported for the fungal enzyme [36]. As 
chlorate was the most potent inhibitor of GmATPSl it was used 
for inhibition studies of the forward reaction. Chlorate was an un- 
competitive inhibitor (K { = 132 ± 4 /xM) versus ATP (Figure 3 A) 
and competitively inhibited (K { = 103 ± 19 /xM) GmATPSl with 
respect to sulfate (Figure 3B). The observed inhibition patterns 
indicate that ATP binds first, followed by sulfate in the forward 
reaction [38]. 



Reverse ATP synthesis reaction: initial velocity 
analysis 

Using the bi substrate variation approach described for the for- 
ward reaction, the ATP synthesis or reverse reaction catalysed 
by GmATPS 1 was examined. Initial velocity patterns for the en- 
ergetically favourable reverse reaction showed intersecting lines 
consistent with a sequential type of mechanism (Figure 4). Global 
fitting of the bi substrate data for the reverse reaction yielded a 
best fit to an ordered sequential mechanism (r 2 = 0.961) com- 
pared with a random sequential mechanism (r 2 = 0.909). The 
fitted kinetic parameters for the reverse reaction are summarized 
in Table 1 and are similar to previously reported values [10]. 

Reverse ATP synthesis reaction: ITC analysis 

To determine the order of binding, and thus the kinetic mech- 
anism, in the reverse reaction, ITC was used. Product inhibi- 



tion by ATP was complicated by the use of a coupling system 
and the binding of oxyanions (sulfate, chlorate, phosphate and 
phosphonoformate) was too weak to compete with the energet- 
ically favourable ATP synthesis reaction. Ligand-binding ana- 
lysis provides a direct method for distinguishing between random 
versus sequential order of addition. 

Titrations of GmATPS 1 with either APS or PPi yielded clearly 
distinct results (Figures 5A,B). Addition of APS showed an exo- 
thermic heat signature indicative of tight binding to GmATPS 1 . 
Fitting of the data to a one- site binding model provides a de- 
termination of the K d value that is comparable with the K m of the 
reaction (Table 2). Titrations using PPi as the ligand did not indic- 
ate interaction with GmATPSl, as no heat signature was detected 
(Figure 5B). These results were consistent with initial binding of 
APS, followed by PPi. To confirm this order of events, GmATPS 1 
was pre-incubated with APS, and then the E»APS complex was 
titrated with PPi. To prevent catalysis, Mg 2+ was not added to 
either the ligands or protein [36]. As shown in Figure 5(C), this 
experiment resulted in a clear heat signature and allowed for the 
determination of thermodynamic parameters for formation of 
the ternary complex (Table 2). This suggests that formation 
of the E»APS complex is necessary for PPi binding. 



DISCUSSION 



Sulfur, like nitrogen, phosphorus and potassium, is a macronutri- 
ent essential for plant growth, crop yields and resistance to patho- 
gens [1-5,39]. In plants, ATP sulfurylase catalyses the thermo- 
dynamically unfavourable formation of a mixed phospho sulfate 
anhydride in APS as the first committed step in assimilation of 
environmental sulfur [8-1 1]. As part of this chloroplast-localized 
pathway, the activity of ATP sulfurylase needs to be coordinated 
with downstream enzymes to control flux through the system for 
maintaining sulfur supplies and to ensure that its products (APS 
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Figure 5 ITC analysis of ligand binding to GmATPSI 

(A) Representative experimental data for the APS titration is plotted as 
heat signal (/xcal s _1 ) versus time (min). Each experiment consisted 
of 20-30 injections of 10 fi\ each of ligand into a solution containing 
protein. (B) Titration of GmATPSI with APS (solid squares) and PPj 
(open circles). The integrated heat response is shown with the solid 
line representing the fit of the data to a one-site binding model. (C) 
Titration of the GmATPSl»APS complex with PPj. The integrated heat 
response is shown with the solid line representing the fit of the data to 
a one-site binding model. 



and PPi) do not accumulate [10]. In contrast with the multi- 
functional, allosterically regulated ATP sulfurylases from mam- 
mals, fungi and bacteria, the plant enzyme functions as a dimeric 
protein and has not been as well studied biochemically [10,24- 
28]. Since engineering of plant sulfur assimilation has a number 
of potential biotechnology applications, understanding the bio- 
chemical control and chemistry of the enzymes in this pathway is 
useful for manipulating the pathway [4-6,39-42]. To better un- 
derstand the molecular function of the plant ATP sulfurylase, we 
used a combination of initial velocity, inhibition and ITC-binding 
experiments to examine the kinetic mechanism of this enzyme. 

Physiologically, ATP sulfurylase catalyses the formation of 
APS and PPi from ATP and sulfate (Figure 1). Kinetic analysis 



of the forward reaction of GmATPSI indicates that a single- 
displacement or ternary complex mechanism occurs (Figure 2). 
The intersecting families of lines for each substrate eliminates 
consideration of a ping-pong type of reaction, but simultaneous 
global data fitting did not adequately distinguish between the 
sequential versus random order kinetic mechanisms. For the in- 
hibition studies of the forward reaction, the use of chlorate as a 
dead-end inhibitor was chosen over product inhibition analysis 
with APS and PPi because of how tightly the products inhibit the 
enzyme and the effect of APS on the coupling enzymes in the as- 
say [8,36]. The competitive inhibition versus sulfate and uncom- 
petitive inhibition versus ATP observed for chlorate (Figure 3) 
is consistent with ATP binding first in an ordered mechanism. 
Based on the dead-end inhibition data, the kinetic parameters for 
an ordered sequential kinetic mechanism are reported in Table 1 . 

For the reverse or ATP synthesis reaction catalysed by 
GmATPSI, the bi substrate variation experiment agrees with the 
single-displacement mechanism observed in the forward reac- 
tion (Figure 4). In this case, global data fitting yielded a better 
agreement with an ordered sequential mechanism compared with 
a random ordered mechanism (Table 1). To distinguish between 
the two possible kinetic models, ITC was used to directly ex- 
amine APS and PPi binding to GmATPSI (Figure 5). This ap- 
proach could not be used in the forward reaction direction because 
of the lack of heat signature associated with binding of either 
ATP or sulfate to GmATPSI (results not shown). Formation of 
GmATPSI binary complex with APS showed tight-binding af- 
finity and an exothermic interaction (Figure 5B and Table 2). In 
contrast, no heat signature was detected in titrations using PPi 
(Figure 5B and Table 2), even though both substrates for the re- 
verse reaction have comparable kinetic parameters (Table 1 and 
[10]). To confirm binding of PPi to the E»APS complex, titration 
to the complex was performed in the absence of magnesium to 
prevent catalysis [36]. Unlike the PPi titration to free enzyme, 
a strong heat signature was detected for formation of the bin- 
ary complex (Figure 5C and Table 2). These results suggest an 
ordered sequential mechanism for the reverse reaction. 

An overall kinetic mechanism for GmATPSI, and the related 
plant enzymes, is now proposed (Figure 6). In the forward re- 
action, an ordered mechanism in which sulfate binding follows 
initial binding of ATP is most consistent with the initial velocity 
and dead-end inhibition data. For the reverse reaction, the ini- 
tial velocity analysis and ITC-binding studies indicate that APS 
binding occurs first, followed by addition of PPi. 

To date, kinetic studies of the ATP sulfurylases from mammals, 
fungi and plants have led to two different models for binding of 
substrates [8,36,43]. Initial velocity and inhibition studies of the 
enzyme from rat and P. chrysogenum were consistent with an 



Table 2 Thermodynamic parameters of ligand binding to GmATPSI 

All titrations were performed at 25 °C as described in the Experimental procedures. ITC data were fit to a one-site binding 


model. ND = not detected. 


Protein in cell 


Titrant 


N 


K d (fiM) 


AG (kcal mol 1 ) 


AH (kcal mol 1 ) 


-7 AS (kcal mol 1 ) 


GmATPSI 
GmATPSI 
GmATPSl.APS 


APS 

PPi 

PPj 


1.03 ±0.01 
ND 

0.96 ±0.01 


0.44 ±0.02 
ND 

0.30 ±0.02 


-8.66 ±0.42 
ND 

-8.89 ±0.51 


-14.4 ±0.1 
ND 

-15.8 + 0.1 


5.69 

ND 

6.92 
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Figure 6 Model for the ordered sequential kinetic mechanism of 
plant ATP sulfurylase 



ordered addition of substrates in each reaction direction [36,43]. 
In contrast, studies with the spinach enzyme suggest a random 
addition of substrates in the forward reaction and an ordered 
addition in the reverse reaction [8]. The kinetic analysis of the 
soybean enzyme presented here and earlier work with the spin- 
ach enzyme [8] differ in the proposed models for the forward 
reaction; however, analysis of the spinach enzyme used a molyb- 
dolysis reaction and not the physiologically relevant chemistry. 
Moreover, it is possible that the random mechanism exhibits a 
preferential binding order, which appears as an ordered mechan- 
ism using steady-state kinetic approaches. For example, with the 
plant APS kinase, steady-state kinetic studies suggested that 
the enzyme used an ordered mechanism for catalysis [44]. Later, 
extensive analysis of ligand binding by ITC demonstrated that 
each substrate interacts with the free enzyme to form binary 
complexes (i.e. random addition), but that the order of substrate 
addition dramatically alters affinity for the second ligand to form 
the ternary complex [30-32]. As noted above, global fitting 
of the initial velocity data for the forward reaction catalysed by 
GmATPSl did not distinguish between the two possible sequen- 
tial mechanisms and the interaction factor obtained for a random 
mechanism suggested that binding of one substrate decreases 
affinity for the second. Nevertheless, the dead-end inhibition pat- 
tern observed with GmATPSl using chlorate supports an ordered 
sequential mechanism. 
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